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Abstract 

We study the production of W bosons at large transverse momentum in pp collisions. 
We show that the next-to-leading order cross section at large transverse momentum is 
dominated by threshold soft-gluon corrections. We add next-to-next-to- leading-order soft- 
gluon corrections to the exact next-to-leading-order differential cross sections. We find 
that these higher-order corrections provide modest enhancements to the transverse mo- 
mentum distribution of the W at the Tevatron, and reduce significantly the dependence 
on the factorization and renormalization scales. 



1 Introduction 



The production of W bosons in hadron colliders is a process of relevance in testing the Standard 
Model and in estimates of backgrounds to new physics. For example, Wbb production is 
the principal background to the associated Higgs boson production pp H{-^ hh)W at the 
Tevatron \X.. 

The calculation of the complete next-to-leading-order (NLO) cross section for W hadropro- 
duction at large transverse momentum was presented in Refs. |21E1- The NLO results in [2||2] 
showed an enhancement of the differential distributions in transverse momentum Qt of the W 
boson. The Qt distribution falls rapidly with increasing Qt, spanning five orders of magnitude 
in the 30 GeV < Qt < 190 GeV region at the Tevatron. The NLO corrections also considerably 
stabilize the dependence of the cross section on the factorization and renormalization scales. 

The calculation of hard-scattering cross sections near partonic threshold, such as electroweak- 
boson production at high transverse momentum, involves corrections from the emission of soft 
gluons from the partons in the process. At each order in perturbation theory one encoun- 
ters large logarithms that arise from incomplete cancellations near partonic threshold between 
graphs with real emission and virtual graphs. This is due to the limited phase space available 
for real gluon emission near partonic threshold. These threshold corrections, calculated in the 
eikonal approximation, exponentiate as a result of the factorization properties |31 El E] of the 
cross section. The cross section is factorized into functions that describe gluons collinear to 
the incoming partons, hard quanta, and noncollinear soft-gluons. The renormalization group 
properties of these functions result in resummation. Such threshold corrections have by now 
been successfully resummed for many processes |7j. 

In this paper we further the study of the resummation of threshold logarithms and the ex- 
pansion of the resummed cross section at next-to-next-to-leading order (NNLO) for electroweak 
boson hadropro duct ion, which was first presented in Ref. at next-to-next-to-leading loga- 
rithmic (NNLL) accuracy. In Ref. no numerical phenomenological studies were made. In 
this paper we study the significance of NNLO soft-gluon corrections for ly-boson production at 
large transverse momentum. Furthermore we increase the accuracy of the theoretical calcula- 
tion to next-to-next-to-next-to- leading logarithms (NNNLL). Related studies for direct photon 
production, including numerical results, were made in Refs. jHlEI. The partonic subprocesses 
involved in direct photon production are similar to the ones discussed in this paper. 

Recently, a unified approach to the calculation of NNLO soft and virtual corrections for 
processes in hadron-hadron and lepton-hadron colliders has been presented in Ref. fl]. This 
work unifies and extends previous approaches by going beyond NNLL accuracy. We follow 
this reference and calculate the NNLO soft corrections to W production at large transverse 
momentum at the Tevatron. The form of our theoretical results differs slightly from (but is 
consistent with) the one in Ref. |H], partly in that here we use the transverse momentum Qt 
as the hard scale instead of the mass of the W . 

We note that there has also been work on resummation of Sudakov logarithms at small 
transverse momentum for electroweak boson production ^21 El as well as on joint resummation 
of Qt and threshold logarithms . 

In Section 2 we discuss some of the kinematics of the partonic subprocesses involved and 
we define the corrections to be calculated. In Section 3 we provide results for the NLO soft and 
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virtual, and the NNLO soft-gluon corrections. The NLO soft corrections agree near partonic 
threshold with the exact NLO calculations in Refs. [3 El while our NNLO soft-gluon corrections 
provide new predictions. In Section 4 we study the numerical effect of the NNLO soft-gluon 
corrections to W hadroproduction at the Tevatron Run I and II, while also showing that the 
NLO cross section is dominated by soft-gluon corrections. 

2 Kinematics and plus distributions 

For the hadronic production of aW boson, with mass mw, 

JiAiPA) + hBiPn) W{Q) + X , (2.1) 
we can write the factorized single-particle-inclusive cross section as 

Eq ^^fcAfeB--4y(Q)+X ^ ^y"rfxirfX2 0/,//,^(x2,/i^) 



X Eq f^F, (2.2) 



where Eq = Q^, (pf/h is the parton distribution for parton / in hadron h, and a is the pertur- 
bative parton-level cross section. The initial-state collinear singularities are factorized into the 
parton distributions at factorization scale fip, while fin is the renormalization scale. 

At the parton level, the lowest-order subprocesses for the production of a boson and a 
parton are 

q{Pa) + g{Pb) w{Q)+q{p,), 

q{Pa)+qiPb) W{Q)+g{p,). (2.3) 

The hadronic and partonic kinematical invariants in the process are 

S={PA + PBy, T={Pa-Q)\ U={Pb-Q)\ S2 = S + T + U -q\ 

S={Pa+ Pb)\ t={pa- Q)\ U = {pb- Q)^ S2 = S + t + U-Q^, (2.4) 

where 5*2 and S2 are the invariant masses of the system recoiling against the electroweak boson 
at the hadron and parton levels, respectively. S2 = (pa+Pb — Q)'^ parametrizes the inelasticity of 
the parton scattering, taking the value S2 = for one-parton production. Since Xj is the initial- 
parton momentum fraction, defined by pa = xiPa and pb = X2Pb, the hadronic and partonic 
kinematical invariants are related by s = X1X2S, t — Q'^= xi{T — Q"^), and u — Q'^ = X2{U — Q"^). 

In general, the partonic cross section a includes distributions with respect to S2 at n-th 
order in the strong couphng of the type 

"In'"(s2/Q^)l . ^ 
^ , m<2n-l, 2.5 

52 
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defined by their integral with any smooth function / by 



These "plus" distributions are the remnants of cancellations between real and virtual con- 
tributions to the cross section. Note that in Ref. we used Q instead of Qt in the plus 
distributions. Here we prefer to use Qt as we find it a slightly better physical hard scale for 
the Qt distributions that we will be calculating. Below we will make use of the terminology 
that at n-th order in as the leading logarithms (LL) are those with m = 2n — 1 in Eq. ()2.5|) . 
next-to- leading logarithms (NLL) with m = 2n — 2, next-to-next-to-leading logarithms (NNLL) 
with m = 2n — 3, and next-to-next-to-next-to-leading logarithms (NNNLL) with m = 2n — 4. 



f ^2 max 

/ dS2f{s2) 

Jo 



\n"'{s2/Q 



S2 



3 Next-to-next-to-leading order soft-gluon corrections 

In this section, for each of the subprocesses in Eq. (j2.3|) we first present the NLO soft and 
virtual corrections. We then present the NNLO soft-gluon corrections. We work in the MS 
scheme throughout. 

3.1 The qg — > Wq subprocess 

The Born differential cross section for this process is 



da, 



d^Q 



(3.1) 



where 



^qg^Wq 



^qg 



(3.2) 



's t 2uQ^ 
As st 



with L the left-handed couplings of the W boson to the quark line, / the quark flavor and X)/ 
the sum over the flavors allowed by the CKM mixing and by the energy threshold. For the L 
couplings we choose the conventions of Ref. jSj. Also Cp = {N^ — l)/(2A^c) with Nc = 3 the 
number of colors. 

We can write the NLO soft and virtual corrections for qg — > Wq in single-particle inclusive 
kinematics as 



d^Q 



Q- 



pB 

qg^Wq 



TT 



S2 



+ cf 



+ Cf 5{S2) 



(3.3) 
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Note that the [ln(s2/Q|')/s2]+ term (which is the LL, since m = n = 1 in Eq. ()2.5p ) and 
the [1/52]+ term (NLL, since n = 1, m = 0) are the soft gluon corrections. The 6(32) term is 
the contribution from the virtual corrections. Below we use the terminology "NLO-NLL" to 
indicate when, at next-to-leading order, we include the LL and NLL soft-gluon terms (as well 
as scale dependent terms in 6(32)) ■ Also the terminology "soft and virtual" is used to denote 
all the terms in Eq. ()3.3|) . As we will see, we need to know both soft and virtual corrections at 
NLO in order to derive the soft terms at NNLO with at least NNLL accuracy. 

The NLO coefficients in Eq. Q are cf = Cp + 2Ca, 



C2 



(Cp + Ca) In 



'/4 



(3.4) 



and 



^99 
•-1 



„95 

[Br + Bfnj + Cf + Cl^Uf] + ^ InM ^ ) + cf In ( ^ 



'Q 



2Ai9 



'Q 



(3.5) 



with Ca = Nc, Uf = 5 the number of light quark flavors, and Bf^ , i?!^, Cf^ , and Cf^ as 
given in the Appendix of Ref. but without the renormalization counterterms and using 
fA ^ ln(A/Q2) = 0. 

Note that we can write cl^ = C2^ + with C2^ = —{Cf + Ca) ln(/i|./s). Similarly we also 
write Ci = cf^ + Tf ^ with 



^19 

-111 




(3.6) 



where (3q = (IIC^ — 2nj)/3 is the lowest-order beta function. Thus, c^^^ and cf^^ are scale- 
dependent parts of the c^^ and Ci coefficients, respectively, while and Tf^ are scale- 
independent parts. We have kept the factorization scale fip and the renormalization scale 
fj,ji separate. Finally, another useful notation is Ci = c'l^^ + T'f where c'f^ is defined as 
with Q\ instead of 3 in the denominators of the logarithms involving the scales ^p and ^r. 

Using the above conventions, the NNLO soft and virtual corrections for qg — > Wq can be 
written as 



with 



qg^Wq' 



■ a' 



qg^Wq 



(3.7) 



,(2) 

qg^Wq 



ln\32/Q^T) 



S2 



+ Cp 



+ 



.19 „19 



-Gr o 



f^O qg ,r^° 



ln\32/Q^T) 
S2 



2 4 Is 




I J J19 (19 /- „<19 „<19 1 {' 
+ 1^2 Ci — L,2 C2 C3 + (,3 {O 



99 \2 




^ T|5 + ^^,9 In fd\ + (^Cp + Ca) K 



HS2IQI) 



19 
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+ {Cf + Ca) In (^C2Tf - UCf + Ca) In (^^) - Cacf 



32 """ V Q^y 16 



+ 



8 



-C. In 



-4 -y-lD-^ 



17 g/g + 7g/3Jin 



/(2) 



^2 



(3i 



where = C^(67/18 — C2) — 5nj/9 is a two-loop function in the MS scheme [THj, C2 = 7r^/6 
and Cs = 1.2020569 ■ ■ ■ are Riemann zeta functions, (3i = 34C^/3 — 2nf{CF + 5Ca/3) is the 
next-to-leading order beta function, and 



/(2) 
^1/1 



cl 



32 



^ 3 \ ^ ^ / 3^ 11 
C2 + 2C3j + CfCa (^-^Cs + 
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48 



(3.9) 



7'^ 



(3.10) 



are two- loop parton anomalous dimensions jTHl 1^ . 

The soft corrections are the [ln^(s2/Q^)/s2]+ term (which is the LL, since n = 2, m = 3 in 
Eq. (Q), [In2(s2/Q|)/S2]+ term (NIL, n = 2, m = 2), [ln(s2/Q^)/s2]+ term (NNLL, n = 2, 
m = 1), and [l/s2]+ term (NNNLL, n = 2, m = 0). The function ^g) ^^e NNNLL term 
denotes a set of two-loop contributions jTUl E] and is given by 



7(2) 
^qg 



189 



\4 12^ 32 



+ ^1(7C3 + ^C2 



11 



^j+-/^46C2 + Y^ 



:C2 



5 

10^ 



. (3.11) 



Note that we have not included in Q^^^ two-loop process-dependent contributions; however 
from related studies for other processes, including top hadroproduction ^Hl and direct photon 
production 10 we expect such contributions to be small. It is actually the —(2C2^cl^ + C3(c3^)^ 
terms in Eq. (j3.8p that provide the major contribution to the NNNLL term. 

We use the terminology "NNLO-NNLL" below, to indicate that we include the LL, NLL, 
and NNLL terms at NNLO, while we use "NNLO-NNNLL" to indicate that in addition we 
include the NNNLL terms as well. 
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The term proportional to S{s2) includes the virtual corrections. We remind the reader that 
here we do not calculate the full virtual corrections; the term R^J that denotes the scale- 
independent virtual corrections is currently unknown. However, we have calculated explicitly 
here all the 6(32) terms that include the scale-dependence. Note that when discussing scale 
dependence at NNLO-NNNLL we include all scale-dependent 6(32) terms, which is a consistent 
approach from the resummation procedure, see Ref. [T^ . 



3.2 The qq — > Wg subprocess 



Next, we consider the qq — > Wg partonic subprocess. Here the Born differential cross section 
is 



da?. 



Q- 



Fqq^Wg ^(-^s) , 



(3.12) 



where 



qq^Wg 

Am 



u t 20^3 

7 + - + ^ 
t u tu 



j^qq 
2, 



'fbfal , 



(3.13) 



The NLO soft and virtual corrections in single-particle inclusive kinematics can be written 



as 



a CTqq^Wg 

d^Q 



Q- 



^qq^Wg 



TT 



J M 


MS2/QI)' 




■ 1 ■ 


^ + Cf5{32)^ 


r 


■^2 


+ 


-S2- 



(3.14) 



Here the NLO coefficients are = ACp — Ca, 



J11 
^2 



2C^ln(gl -{2Cp-CA)\n 



tu 

3Qj: 



4 



(3.15) 



and 



^11 



Bf + Cf + {Bf + Dl^J] 



n. 



+ ^ln2 
2 



Q't 



+ cf In 



(3.16) 



with B^ ^ B^ ^ 1 and D^^-* as given in the Appendix of Ref. but without the renormal- 
ization counterterms and using = 0. Again, we can write = cf^^ + Tf' with 



nil 



In ^ 



(3.17) 



and = cl''^ + Tl' with c^"^^ = — 2Ci? ln(yu|,/s). Finally, another useful notation is cf' = 
c'f^ + T'f^ where c'f^ is defined as c^^^ with instead of 3 in the denominators of the 
logarithms involving the scales fip and fiR. 
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The NNLO soft and virtual corrections for qq — > Wg can be written as 



da. 



(2) 

qq^Wg 



with 



,(2) 

qq^Wg 



1 



^ qq^Wg ^2 qq^Wg 



(3.18) 



S2 



+ 



+ Ca 



(2 (4 



3 
2 

2 



„m „qq n'^i < r . 



/5o 



ln^(g2/g|) 

S2 



Tf+^cf In 



+ 2Ci.K 



K 
~2 



4 



+ i ^2 Ci — (,2 C2 C3 + (,3 (^Cg 



,99^2 



S2 

In 



16 



+ Ca 



/^O i„2 / 



2 I s 



In 

16 



^ ^qq 



1 

S2 



2CIC2 In^ f ^"j + 2Cf In f ^"j f (2^^" - 2(2^^ In ( ^ 



Qt 



.qq 



^0 1 2 / A'-R 



- i^ln 
32 



^Q2y 16 \Ql^ 



+ 



2 



+ Inl'^-Vn^^^ 



2 i ^^q/q- 



The function denotes again a set of two-loop contributions ITT] and is given by 



,(2) 

^qq 



^ ^ 22^ 299\ ^ / 4^ 50\ 

,(2) 



(3.19) 



(3.20) 



Again, we have not included in Qgg' two-loop pro cess- dependent contributions. 

We also note that we do not calculate the full virtual corrections. All ^(^2) terms shown 
here only include the scale-dependence and are used in the study of the scale dependence of 
the NNLO-NNNLL cross section. The term Rgg is unknown. 



4 Numerical results 



For the W hadroproduction cross section for hA{PA) + hBiPs) W{Q) + X we defined the 
hadronic and partonic kinematical invariants in Section 2. The hadronic kinematical variables 
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T and U can be written alternatively as T = — mxy/Se ^ and U = niy^ — mT\fSe^ . Here 

miy is the W mass, = \J Q\ + is the transverse mass, and y is the rapidity of the vector 
boson. The differential Qt distribution can then be written as 




where F = In (5 + ^^2-1), B = + m^)/(2mTV^), y = Y{2y'-1), 

- VSrriTey + Xi [S - v/S^mTC"^) , (4.2) 

\fSmTey-m'L S2 - mL + \fSmTXie~y ,^ „, 

1= X2 = -1= . (4.3) 

S — \/SniTe~y XiS — V SuiTey 

We now apply our results to W hadroproduction at large transverse momentum at the 
Tevatron. Throughout we use the MRST2002 approximate NNLO parton densities [201 • In 
Fig. 1 we plot the transverse momentum distribution, da/dQ"^, for W hadroproduction at 
the Tevatron Run I with = 1.8 TeV. Here we have set fJ-F = fJ-R = Qt- We plot Born, 
exact NLO 013], NNLO-NNLL, and NNLO-NNNLL results. We see that the NLO corrections 
provide a significant enhancement of the Born cross section. The NNLO-NNLL corrections 
provide a further modest enhancement of the Qt distribution. If we increase the accuracy by 
including the NNNLL contributions, which are negative, then we find that the NNLO-NNNLL 
cross section lies between the NLO and NNLO-NNLL results. Since it is hard to distinguish 
between the curves in Fig. 1, we provide another figure. Fig. 2, which emphasizes the high-Q^ 
region where the soft-gluon approximation holds best. 

In Fig. 3 we plot the K-factors, i.e. the ratios of cross sections at various orders and 
accuracies to the Born cross section, all with lip = I^r = Qt, in the high-Q^ region. We also 
show the ratio of the exact NLO to the NLO-NLL cross section. It is clear from this line being 
very close to 1 that the NLO-NLL result is a very good approximation to the full NLO result, 
i.e. the soft-gluon corrections overwhelmingly dominate the NLO cross section. The difference 
between NLO and NLO-NLL is only 2% for Qt > 90 GeV and less than 10% for lower Qt 
down to 30 GeV. The fact that the soft-gluon corrections dominate the NLO cross section is 
a major justification for studying the NNLO soft gluon corrections to this process. We can 
also see that the various ii"-factors shown in Fig. 3 are moderate, and nearly constant over 
the Qt range shown even though the distributions themselves span two orders of magnitude 
in this range. The NLO corrections are nearly 30% over the Born, with the NNLO corrections 
giving an additional increase, so that the NNLO-NNLL/Born i^'-factor is around 1.4 and the 
NNLO-NNNLL/Born TT-factor is around 1.35. 

In Fig. 4 we plot the scale dependence of the differential cross section for Qt = 80 GeV. 
We define = fip = fiji and plot the differential cross section versus ^i/Qt over two orders 
of magnitude: 0.1 < ^i/Qt < 10. We note the good stabilization of the cross section when 
the NLO corrections are included, and the further improvement when the NNLO-NNNLL 



max 
bo 



A 



9 



pp-> w 



S'"=1 .8 TeV 




30 50 70 90 110 130 150 170 190 

Qj (GeV) 



Figure 1: The differential cross section, da/dQ'^, for W liadroproduction in pp collisions at 
the Tevatron with =1.8 TeV and /iF = /^/? = Qt- Shown are the Born (solid line), NLO 
(long-dashed line), NNLO-NNLL (short-dashed line), and NNLO-NNNLL (dotted line) resuhs. 

corrections (which include all the soft and virtual NNLO scale terms) are added. The NNLO- 
NNNLL result approaches the scale independence expected of a truly physical cross section. 

In Fig. 5 we plot the differential cross section da/dQ'^ at high Qt with y/S = 1.8 TeV for 
two values of scale, Qt/2 and 2Qt, often used to display the uncertainty due to scale variation. 
We note that while the variation of the Born cross section is significant, the variation at NLO 
is much smaller, and at NNLO-NNNLL it is very small. In fact the two NNLO-NNNLL curves 
lie on top of the /i = Qr/2 Born curve, and so does the NLO /i = Qt/2 curve. These results 
are consistent with Fig. 4. 

In Fig. 6 we plot the transverse momentum distribution, da/dQ"^, for W hadropro duct ion 
at the Tevatron Run II with = 1.96 TeV. Again, we have set fip = /^i? = Qt- The relative 
size of the corrections is similar to Fig. 1. In Fig. 7 we display the high-Q^ region more clearly. 
We do not plot the i^-factors for Run II as they are nearly identical to those for = 1.8 
TeV. 

In Fig. 8 we plot the differential cross section da/dQ^ at high Qt with = 1.96 TeV for 
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pp-> w 




I I I I I I I I I I I I I I I I 

100 120 140 160 180 
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Figure 2: The differential cross section, da/dQ'^, of Fig. 1 at liigli Qt- Tlie labels are the same 
as in Fig. 1. 

two values of scale, Qt/2 and 2Qt. The results are analogous to those in Fig. 5. 

5 Conclusion 

We have presented the NNLO soft-gluon corrections for W hadropro duct ion at large transverse 
momentum in pp collisions, in particular at the Tevatron Run I and II. We have shown that 
the NLO soft-gluon corrections completely dominate the NLO differential cross section at large 
transverse momentum, while the NNLO soft-gluon corrections provide modest enhancements 
and further decrease the factorization and renormalization scale dependence of the transverse 
momentum distributions. 
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pp->W K-factors 

S''^=1 .8 TeV 

I 1 1 1 1 1 1 1 1 1 1 1 1 1 

NLO/Born J 

NLO-NLL/Born 

NNLO-NNLL/ Born -- 

NNLO-NNNLL/ Born 
NLO/NLO-NLL 



100 120 140 160 180 

Qj (GeV) 

Figure 3: The i^'-factors for the differential cross section, da/dQ'^, for W hadroproduction in pp 
colhsions at the Tevatron with = 1.8 TeV and /j^f = I^r = Qt- Shown are the K- factors for 
exact NLO/Born (long-dashed line), NLO-NLL/Born (dash-dotted line), NNLO-NNLL/Born 
(short-dashed line), and approximate NNLO-NNNLL/Born (dotted line) results. Also shown 
is the ratio of the exact NLO to the NLO-NLL cross section (solid line). 
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Figure 6: The differential cross section, da/dQ'^, for W hadropro duct ion in pp collisions at the 
Tevatron Run II with = 1.96 TeV and jJ^F = I-J'R = Qt- Shown are the Born (solid line), 
exact NLO (long-dashed line), NNLO-NNLL (short-dashed line), and NNLO-NNNLL (dotted 
line) results. 
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Figure 7: The differential cross section, da/dQ"^, of Fig. 6 at liigli Qt- Tlie labels are the same 
as in Fig. 6. 
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Figure 8: The differential cross section, da/dQ"^, for W hadropro duct ion in pp collisions at the 
Tevatron with = 1.96 TeV and ^ = = = Qt/'^ or 2Qt- Shown are the Born (sohd 
lines), NLO (long-dashed lines), and NNLO-NNNLL (dotted lines) results. The upper lines are 
with /i = Qt/2, the lower lines with ^ = 2Qt- 
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